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Abstract 

A new variational principle - extremizing the fixed frame kinetic 
energy under constant relative enstrophy - for a coupled barotropic 
flow - rotating solid sphere system is introduced with the following 
consequences. In particular, angular momentum is transfered between 
the fluid and the solid sphere through a modelled torque mechanism. 
The fluid's angular momentum is therefore not fixed but only bounded 
by the relative enstrophy, as is required of any model that supports 
super-rotation. 

The main results are: At any rate of spin 0, and relative enstro- 
phy, the unique global energy maximizer for fixed relative enstrophy 
corresponds to solid-body super-rotation; the counter-rotating solid- 
body flow state is a constrained energy minimum provided the relative 
enstrophy is small enough, otherwise, it is a saddle point. 

For all energy below a threshold value which depends on the rel- 
ative enstrophy and solid spin $7, the constrained energy extremals 
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consist of only minimizers and saddles in the form of counter-rotating 
states. Only when the energy exceeds this threshold value can pro- 
rotating states arise as global maximizers. 

Unlike the standard barotropic vorticity model which conserves 
angular momentum of the fluid, the counter-rotating state is rigorously 
shown to be nonlinearly stable only when it is a local constrained 
minima. The global constrained maximizer corresponding to super- 
rotation is always nonlinearly stable. 
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1 Introduction 



The main aim of this paper is to present the simplest possible formulation 
and rigorous analysis of a theory for geophysical flow coupled to an infinitely 
massive rotating sphere, which exhibits characteristics of super and sub- 
rotation. Although the General Circulation Model (GCM) simulations re- 
produced some aspects of super-rotation such as those on Venus and Titan, 
a rigorous theory which can be analysed in classical mathematical terms, is 
lacking. This is the motivation which guides our approach here. We therefore 
propose a mathematically precise constrained variational problem in terms of 
the pseudo energy (rest frame kinetic energy) of a coupled barotropic flow on 
a rotating sphere, in which the only explicit constraints are relative enstrophy 
and zero total circulation. Unlike previous work, the angular momentum of 
the fluid is not flxed. 

As the most promising phenomenological theory to date is one where 
the rotating atmosphere exchanges angular momentum with the planetary 
surface through a complex torque mechanism mediated by Hadley cells, we 
propose a variational theory to find and analyze the nonlinear stability of 
steady-states of the coupled fluid - sphere system. The coupled system is a 
conservative or non-dissipative model in the sense that the energy and angu- 
lar momentum of the combined fluid-sphere system are separately conserved 
in time. We will not go into the details of standard arguments for the validity 
of an inviscid assumption for the interior flow in the context of large scale 
geophysical flows. Suffice it to say at this point that the role of viscosity is 
restricted to very thin Ekman boundary layers, and is within that part of the 
coupled fluid-sphere system modelled by the complex torque mechanism. 

The main results of this approach is that, at any rate of spin Q and rel- 
ative enstrophy Qrei, the constrained global energy maximum for fixed rela- 
tive enstrophy corresponds to solid-body rotation, w^Ma.,:{Qrei) = VQreii^io, 
in the direction of fl. Another solution, the counter-rotating steady state 
''^miniQrei) — — VQre^V'iO' ^ Constrained global energy minimum provided 
the relative enstrophy is small enough, that is, Qrei < ^^C^ where C = 1 1 cos 6 
II2. If r^^C^ < Qrei < 4ri^C^, then w'^^^{Qrei) is a saddle point. And if 
Qrei > 4Q^C^, then w^^^^{Qrei) is again a saddle point, but in the restricted 
sense that it is a constrained local energy maximum in all eigen-directions 



3 



except 'ipi ±1, '^miniQrei) IS a constramed local energy minimum in the projec- 
tion onto span{ipioj i'l ±i) where the ipi ±i eigenmodes are relative vorticity 
patterns associated with the tilt instability. 

The main result can be viewed another way. For all energy H below a 
threshold value He which depends on the relative enstrophy Qrei and spin Q, 
the constrained energy extremals consist of only minima and saddles in the 
form of counter-rotating states, —y/Qreii^w Only when the energy exceeds 
this threshold value He, can pro-rotating states, VQrei'^io arise as constrained 
global energy maxima. 

It is also useful to view the approach in this paper in the light of a verti- 
cally coarse-grained variational analysis of the steady states of the barotropic 
component of fully stratified rotating flows coupled to a massive rotating 
sphere. We recall another robust connection between the time-independent 
variational theory presented and dynamic initial value problems for damped 
driven quasi-geostrophic flows: The Principle of Selective Decay or Minimum 
enstrophy states that, for suitable initial data, the long time asymptotics of 
the 2D Navier-Stokes equations - damped and forced on a doubly-periodic 
domain (and by extension the sphere) - is characterized by a decreasing 
enstrophy to energy ratio |6j. These Minimum Enstrophy flows are charac- 
terized by special steady states which are the solid body rotational flows in 
the case of the sphere, by virtue of an application of Poincare's Inequality 
j2S|- In 1953, Fjortoft jS] used a spherical harmonics expansion to show that, 
for dissipative barotropic flows, because the enstrophy spectrum is related to 
the energy spectrum by Z{l,t) = l{l + l)E{l,t), the energy fluxes towards 
low wavenumbers / while enstrophy fluxes towards high / modes. His work 
is influential in bringing forth the concepts of the inverse cascade of energy 
and forward cascade of enstrophy that are closely related to Selective Decay. 

Indeed, an inviscid asymptotic formulation of the Principle of Minimum 
Enstrophy yields a variational problem that is the dual of the one in this pa- 
per, namely, extremize enstrophy at fixed values of the total kinetic energy 
of flow without fixing the angular momentum of the fluid. The correspond- 
ing extremals including states of minimum enstrophy are again solid-body 
rotating flows. Leith formulated a related method precisely by looking 
for constrained enstrophy extremals in vortex dynamics with fixed angular 
momentum of the fluid. 

Stability results for barotropic flows are discussed in Fjortoft [7j, Tung 
[221, Shepherd et al [21, [ID] amongst many others. Fjortoft [7j first used 
integral theorems to argue that solid-body rotational fiows are always stable 
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in the context of the Hamiltonian PDE known as the barotropic vorticity 
equaton (BVE). The discrepancy between his result and those reported here 

- the retrograde sohd-body flow is always a stable steady-state in the BVE 
but has only conditional nonlinear stability in the coupled barotropic flow 

- rotating solid sphere system - is partly due to the simple but vital fact 
that the fluid's angular momentum is allowed to vary in the latter. We will 
elaborate on this point below. 

Arnold jTTj, ^2] formalized this body of results known as the energy- 
Casimir method, and proved that they work at finite amplitudes, that is 
in a nonlinear way. According to an argument of Shepherd in his refined 
application of the energy-Casimir method to the BVE the Lyapunov 
function method fails precisely when the planetary spin Q is small - the 
proof of the nonlinear stability of solid-body flow states in the BVE breaks 
down when the spin is small or zero. See also the note by Wirosoetisno 
and Shepherd PU] which proved that there are two important cases in 2D 
Euler flows where the Arnold stability method fails, and one of them is 
the non-rotating sphere. Our result on the conditional instability of the 
retrograde (east to west) solid-body flow and the nonlinear stability of the 
prograde sohd-body flow is partially compatible with Shepherd et al because 
we proved instability of the retrograde flow state only when the planetary 
spin fl is smaller than ^jQrei/C^, and because our result also states that the 
prograde rigid rotation state is nonlinearly stable even when the planetary 
spin is relatively small. Again the discrepancy between our results and those 
of Shepherd et al can be traced to the fact that the coupled barotropic fluid 
sphere model is not - unlike the BVE - a Hamiltonian PDE and angular 
momentum of the fluid is not conserved. 

Indeed, the correct variational theory proposed here is related to the 
zero temperature version of the statistical mechanics energy - relative en- 
strophy theory given in [2], ^Hl- There as here, the underlying model - a 
coupled fluid sphere system - cannot be represented by a Hamiltonian or a 
Lagrangian. Therefore, there is no local in time equation of motion for the 
underlying model but rather a generalized variational principle operating 
in overall phase space determines its symmetries and dynamical properties. 
This variational theory and corresponding statistical mechanics theory com- 
prise a new example of Feynman's generalized Least Action Principle that 
is applicable to non-local situations which do not have a Hamiltonian or La- 
grangian formulation. In this setting, we propose a non-Hamiltonian energy 
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functional that plays a principal role in the theory presented here and, as the 
integrand in the action of a partition function (or path- integral) 0, in the 
corresponding statistical mechanics theory as well. 

This statistical mechanics theory predicts that the macrostates with ex- 
tremal internal energy - highly organized, maximally solid-body flow states 
- are invariably associated with low Shannon entropy, and are therefore the 
most probable states when the absolute value of the statistical temperature is 
small. On the other hand, when the temperatures have large absolute values, 
the most probable macrostates - disordered flow states with little or zero net 
angular momentum - have intermediate values of internal energy and high 
Shannon entropy 0. 

In section 2 we give properties of the model, as useful background for the 
rest of the paper. Section 3 describes the energy-enstrophy class of variational 
formulations for the coupled fluid - sphere model. Section 4 gives the con- 
strained optimization of flow energy, H on fixed relative enstrophy manifolds 
and section 5 describes the physical characteristics of the constrained energy 
extremals that arise in the model. Section 6 gives the nonlinear stability 
results. Section 7 concerns physical consequences and concluding remarks. 

2 Coupled Barotropic Fluid - Rotating solid 
Sphere Model. 

Consider the system consisting of a rotating massive rigid sphere of radius R, 
enveloped by a thin shell of non-divergent barotropic fluid. The barotropic 
flow is assumed to be inviscid, apart from an ability to exchange angular mo- 
mentum and kinetic energy with the infinitely massive solid sphere through 
a complex torque mechanism. We also assume that the fluid is in radiation 
balance and there is no net energy gain or loss from insolation. This provides 
a crude model of the complex planet - atmosphere interactions, including the 
enigmatic torque mechanism responsible for the phenomenon of atmospheric 
super-rotation - one of the main applications motivating this work. 

For a geophysical flow problem concerning super-rotation on a spherical 
surface there is little doubt that one of the key parameters is angular mo- 
mentum of the fluid. In principle, the total angular momentum of the fluid 
and solid sphere is a conserved quantity but by taking the sphere to have 
infinite mass, the active part of the model is just the fluid which relaxes by 
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exchanging angular momentum with an infinite reservoir. It is also clear that 
a 2d geophysical relaxation problem such as this one will involve energy and 
enstrophy. The rest frame energy of the fiuid and sphere is conserved. Since 
we have assumed the mass of the solid sphere to be infinite, we need only 
keep track of the kinetic energy of the barotropic fiuid - in the non-divergent 
case, there is no gravitational potential energy in the fiuid because it has 
uniform thickness and density, and its upper surface is a rigid lid. 

In a nutshell, we need to find a suitable set of constraints for the obvious 
choice of objective functional, namely rest frame kinetic energy of fiow in the 
coupled model. The choice of auxiliary conditions or constraints is not apriori 
obvious, as different choices can have slightly different physical consequences 



We will use spherical coordinates - cos^ where 6 is the colatitude and 
longitude 0. The total vorticity is given by 



where 2f2cos^ is the planetary vorticity due to spin rate Q and w = Aip is 
the relative vorticity given in terms of a relative velocity stream function i/j 
and A is the negative of the Laplace-Beltrami operator on the unit sphere 



2.1 Inverse of the Laplacian 

Key to the analysis in this paper is the inverse integral operator defined by 



m 



q(t; cos9, (f)) = Atp + 2Q cos 9 



(1) 




on S"^. In terms of G, the solution of 



is given by 



A-^(/) 



G[f]{x) 
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2.2 Eigenfunctions of G 



The variational analysis of this problem is based on the fact that there is 
an orthonormal basis for L^(5'^) of eigenfunctions of the Laplace-Beltrami 
operator. And G being the inverse of A, its eigenfunctions consist of the 
spherical harmonics, 

i^lm^ ^ = -I, -J (2) 

with eigenvalue A;^ = + 
Thus, a relative vorticity field, by Stokes theorem, has expansion 

W{x) = «/m^im(^)- (3) 

l>l,m 

A key property that will be established later is that the mode aioip^Q{x) 
contains all the angular momentum in the relative flow with respect to the 
frame rotating at the fixed angular velocity of the sphere. 



2.3 Physical quantities of the coupled barotropic vor- 
ticity model 

The rest frame kinetic energy of the fluid expressed in a frame that is rotating 
at the angular velocity of the sohd sphere is 



HtIq] 



dx 



S2 



dx 



S2 



+ 



2 J52 



dx 



— / dx ibq -\ — / dx ui 



where the zonal and meridional components of the relative ve- 

locity. Up is the zonal component of the planetary velocity (the meridional 
component being zero since planetary vorticity is zonal), and -0 is the stream 
function for the relative velocity. 

Since the second term ^ Jg2dx is fixed for a given spin rate ^2, it is 
convenient to work with the pseudo-energy as the energy functional for the 
model, 

H[w] — — - I dx ipq = — f dx ip{x) [w{x) + 2f2cos6'] 
2 J s'^ 2 J 
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= — - J dx iIj{x)w{x) — Q J dx ip{x) cos6 

= ~{w,G[w])-nC{ij,,,G[w]) 

where C = || cos^||2. Later, it will be convenient to perform one more minor 
adjustment of the energy functional in Lemma 3 to turn it into a positive 
definite quadratic form in the Fourier coefficients of the spectral expansion 
of the relative vorticity w in the spherical harmonics ipijn- The following 
calculation using the self-adjoint and eigenvalue properties of the integral 
operator G, 

-C{^,o,G[w]) = ~C{G[^,,],w) 
C 

1 



shows that 



, dx wix) cosOix) 



A[w]^-C{^,„G[w])- 



2 

is the variable net angular momentum of the flow associated with the relative 
vorticity w{x) (assuming unit mass density). Since 

1 If 

E[w] = — - {w,G[w]) = ~2 J '^^ '4'{x)w{x) 

is the kinetic energy of relative motion, this suggests that the pseudo-energy 
H[q] = — I Jg2 dx ipq has the form 

H[q] = E[w] + nA[w] 

of an energy-momentum functional used in standard variational methods for 
rotating problems. This is however only a coincidence. 

A key difference between our approach and the standard variational anal- 
ysis of rotating problems is that neither angular momentum A[w] nor kinetic 
energy H is fixed in our approach - the coefficient aio in the expansion Q 
is not fixed. 

Total circulation in the model is fixed to be J q dx = J wdx = 0, which 
is a direct consequence of Stokes theorem on a sphere and has nothing to do 
with the fact that the barotropic fiow is assumed to be inviscid. It is easy 
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to see that the kinetic energy functional H is not well-defined without the 
further requirement of a constraint on the size of its argument, the relative 
vorticity field w{x). A natural constraint for this quantity is therefore its 
square norm or relative enstrophy so as to carry out the variational analysis 
in the Hilbert space, L^(5'^). But before we settle on this constraint for the 
variational theory, we consider the total enstrophy of the fiow, 

r[gl = / dx q^. 

Expanding (cf. [T]), 

r\q] = I dx = I dx\w + 2ncos9? 

= I dxw^ + An [ dxwcose + 4^]2C^ (4) 

we find that the last term is a conserved quantity because it is the square of 
the norm of the fixed planetary vorticity, the second term is proportional 
to the variable net angular momentum of the fiuid (relative to the rotating 
frame), and the first term is the relative enstrophy. 

From (jH), we deduce that, unlike the BVE, at most one of the two quanti- 
ties - total enstrophy and relative enstrophy - can be conserved in the coupled 
fiuid-sphere model. In other words, if the total enstrophy is conserved then 
relative enstrophy is not because net angular momentum is not conserved; 
if on the other hand, relative enstrophy is conserved, then total enstrophy 
must change with net angular momentum. In fact it is not possible - from 
what is assumed in the model - to exclude the third possibility that neither 
the relative entrophy nor total enstrophy is conserved. 

Thus, we see that our choice of the relative enstrophy constraint, although 
it is natural and required for a rigorous variational analysis, is not a conse- 
quence of its invariance. Instead the argument that to an extent justifies 
this choice is based on the Principle of Minimum Enstrophy. This Principle 
states that only the ratio of enstrophy to energy - not their separate values - 
should be relevant to the analysis of a quasi-steady state in 2d fiows and this 
ratio should have the minimum allowed value for the given geometry. This is 
clearly dual and equivalent to characterizing quasi-steady states in 2d fiows 
in terms of their energy on iso-enstrophy manifolds. 

The second term in is equal to 4f2 times the variable angular momen- 
tum density of the relative fiuid motion and has units of m'^/ s. The physical 
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angular momentum, given by 



p dx w cos6 = p {w,cos9) , (5) 
J 52 

implies that the only mode in the eigenfunction expansion of w that con- 
tributes to its net angular momentum is aioV'io where ip^Q = a cos 9 is the 
first nontrivial spherical harmonic; it has the form of solid-body rotation 
vorticity. 

All other moments of the vorticity / dx g" are deemed less important in 
many physical theories of fiuid motions including the absolute equilibrium 
statistical mechanics model ^S] and the variational problems below [T7] . 
They are, however, not totally irrelevant |24j . 



3 Energy-relative enstrophy variational the- 
ory 

One of the aim of this approach is to answer the question: For which relative 
vorticity field w{x) is the fiuid's rest frame kinetic energy extremal under 
the single additional constraint that the square- norm of w{x) or relative 
enstrophy is fixed? It turns out to be physically relevant that the Wo{x) that 
extremizes H[q\ also maximizes the net angular momentum A[w] since then, 
these extremal steady states are super or sub-rotating according to the sign 
of aio in the expansion Q of wo{x). 

A similar analysis for the non-rotating case {Q = 0) shows below that 
H[w] = E[w] is maximized under fixed enstrophy by pro-rotating solid-body 
flow wq{x) = al()^plQ{x) for aio > which at the same time maximizes net 
angular momentum A[wo]. This implies immediately the proposition: 

Proposition: For fixed Q > 0, H[q\ is maximized under fixed relative 
enstrophy by the solid - body flow wo{x) which also maximizes relative or net 
angular momentum A[wo]- 

What is not obvious and therefore requires the substantial analysis pre- 
sented in this paper is the asymmetry between the energy maximizer and mn- 
imizer - the counter-rotating solid-body flow wi{x) = awipioix) for aio < 0, 
corresponding to minimal angular momentum, and maximizing E[w] (since 
the energy E[w] is even in w), is a. minimizer of H[q\ for values of relative 
enstrophy that are small compared to the rotation rate Q, but becomes a 
saddle point for larger values of relative enstrophy. 
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Not only are the higher vorticity moments / dx g" for n > 2 ignored in 
this paper (cf. also We also do not constrain the total enstrophy 

dx + A^l dx w cos 6 

nor the angular momentum 

dx w cos 9. (6) 



/52 

The only constraint in this paper besides the zero total circulation condition, 

dx w = 0, (7) 

52 

is that the relative enstrophy, 

/ dx w"^ = Qrel > 0, (8) 

752 

is fixed. 



3.1 Constraint Vrei 

We first analyze the physical consequences of the constraint ||w||2 = Qrei > 
in the definition of the subset 

Vrel = S^We L\S^) I Trel[w] = \\w\\2 = Qrel > 0, TC = J^^ dx W = oj, 

which leave the angular momentum (jU)) constrained only by an inequality. 
Fixing the relative enstrophy and relative circulation, this model allows in 
principle energy extrema which can have up to a maximum amount of super 
or sub-rotation. It has the interesting consequence that angular momentum 
is effectively constrained by an inequality. 

Lemma 1: Let the circulation and relative enstrophy of w ^ L'^{S'^) be 
constrained by ^ and ^ respectively. Then, the angular momentum ^ is 
constrained by an inequality, that is, 

< (9) 

where C > is a constant that does not depend on w. 



dx w cos 9 



52 



12 



Proof: By the Cauchy-Schwartz inequality [231, the result is proved with 



It is not surprising that solid-body rotation vorticity in the form of spher- 
ical harmonic ipiQ = bcos6 is the only zero circulation relative vorticity field 
that maximizes angular momentum for given enstrophy. This is stated pre- 
cisely below: 

Lemma 2: The upper bound on the angular momentum in (0j is achieved 
only for relative vorticity field 




(10) 



QED. 



w = k cos 6 



where k is a constant. 



Proof: When w = kcosO, 




dx w cos 6 



k 




and 




Thus, 




dxWCOsO =CjQrel- 




Conversely, if 



w = k cos 6 + atpijn 



where Im 10 and a 7^ 0, then 




dx wcosO < C\\w 



QED. 
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3.2 Augmented energy functionals 

The Lagrange Multiplier method for Hilbert space gives necessary con- 
ditions for constrained extremals in the form of Euler-Lagrange equations. 
More importantly, the explicit form of the Lagrange Multipliers in terms of 
spin, energy and relative enstrophy, provides the detailed physical relation- 
ships that we seek in order to answer the questions raised. We also extend the 
Lagrange Multiplier method to give a geometrical proof of the existence and 
nonlinear stability of the constrained energy extremals where the constraint 
is fixed relative enstrophy. 

The unconstrained optimization of the augmented energy-relative enstro- 
phy functional arising from the Lagrange Multiplier method yield necessary 
conditions for extremals of the constrained optimization problems. Necessary 
conditions for the extremals of the augmented functionals are in turn given in 
terms of their Euler-Lagrange equations (or Gateaux derivative) which take 
the form of interesting inhomogeneous linear equations involving the inverse 
G of the Laplace-Beltrami operator. Bifurcation values of the multipliers can 
be read off the spectrum of G. Extremals of the augmented energy functional 
typically change type for different regimes of the Lagrange multiplier values 
which are marked by bifurcation values. 

The relative enstrophy constraint gives rise to the following constrained 
optimization problems: 

extremize H[w] on (11) 

Vrel = f^We L\S^) I Trel[w] = \ \w\\2 = Qrel > 0, J^JxW = o| .(12) 

By the method of Lagrange Multipliers (see [21], 1201); the energy-relative 
enstrophy functional 

Erel[w; n, R] = H[W] + XTrel[w] (13) 

is the augmented objective functional for the unconstrained optimization 
problem coresponding to the above constrained optimization problem. 

3.3 Lagrange Multipliers 

The Euler-Lagrange Multiplier theorem states: Let wq G Vrei be an extremal 
of H[w]. Then at least one of the following holds: 
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(1) 6Trel{Wo) = 0, 



(2) 6EreliWo) = SH{Wo) + XSTreliWo) = 0. 

The right way to use this theorem is: (a) find a set of relative vorticity 
w G Vrei which satisfies 6Trei{w) = 0, (b) find a second set of relative vorticity 
w G Vrei which satisfies 5H{w) + XdV^eiiw) = 0, for some A, (c) the set of 
constrained extremals of H[w] is contained in the union of the sets in (a) and 
(b), and (d) the value of A is determined from the value of the fixed constant 

'^rellw] = Qrel > 0. 

In implementing this procedure we first compute the variation or Gateaux 
derivative 

5Trei(w, Aw) = 2 wAw dx. 

By choosing Aw = w, we show that this variation never vanishes for any w, 
i.e., 

(1) 5Treliw)^0. 

The Euler-Lagrange Multiplier theorem then states that any extremal Wq of 
the constrained variational problem (jllj) and (jl2j) must satisfy 

(2) 5EreliWo) = 6H{Wo) + X6Trel{Wo) = 

for some value of A. 

The vanishing of the Gateaux derivative of the augmented functional Erei 
gives the Euler-Lagrange equation for the unconstrained problem, which is 
solved in the next section. 

4 Extremals of the augmented energy func- 
tional 

The analysis below will be based on the spherical harmonics (j2I) which are 
eigenf unctions of the Laplace-Beltrami operator on the sphere. In particular, 
the eigenfunction t/'^q = a cos 6' which has eigenvalue Aio = —2, will play a 
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special role in the characterization of the zonal steady states. Using the 
linearity of G, and the eigenf unction expansion 



(^Imi^lm, (14) 

l>l,m 

of relative vorticity w satisfying Jg2 dx w — 0, the Hamiltonian functional is 
expanded in terms of the orthonormal spherical harmonics, to yield 

H = ~{w,G[w])-nC{^,„G[w]) 

- -It. Y^ + l^Ca^o. (15) 

^ l>l,m ^'"^ ^ 

It follows that the coupling between the relative motion and the planetary 
vorticity occurs through the cigenfunction t/j^Q, which has the form of a vor- 
ticity pattern corresponding to solid-body rotational flow. 

For obvious technical reasons, and without loss of generality, wc will use 
the following lemma to change the Hamiltonian by a constant to a quadratic 
form. 

Lemma 3: For fixed spin Q > 0, the energy H for w satisfying 




is modulo the constant 

H, = -\n'C^ (16) 

equal to the positive definite quadratic form (which we will again denote by 
H) 

1 a? 
- ^ E Y^- (17) 

Proof: By completing the square, 

^ l>l,m "^li^ ^ 



-H" = -7 [aio 



+ A 



ail + ^ 



1-1 
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Im 



1 2 Ir^^ 1 >^ a 

-a^o + -nCaio- - 2^ . 

^ ^ ^ l>l m ^'»^ 



+ 



1 r 



2 2 



E 



2 Z>l,m -^'"^ 



2 2 



2 -, A/to 4 



QED. 

Since all the extremals below have the form kip^^ of solid-body rota- 
tion, it is useful to state the following result. The simple proof is left to the 
reader. Graphs of (fTS|) and (fT^ are depicted in figures 1 and 2. 

Lemma 4: The energy and relative enstrophy of the extremals = 
aioV'io tO'kes the form 



^^[aioV'io] ■ 

Qrel = ^rel[oilo4'lO 

Furthermore, for fixed H, 



(18) 



a 



10- 



with solutions 



Qrel ± 2nC^Qrel + ^^C^ = 4:H 

Qrel = (±nC + ^/IH)\ 



(19) 
(20) 



4.1 Lagrange Multipliers 

This constrained variational model is based on the fixed relative enstrophy 
constraint Vrei- By a theorem in [21] on the Lagrange Multiplier method, the 
variational problem (jlip . (jl2j) can be reformulated in terms of extremals of 
the augmented functional, 

Erel[w;Q] = H[w] + XrelTrel[w], 
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Expanding Fre/fw] in terms of (|14j) and combining with (jl7p yields 



a,. 



^ l>l,m ^l"^ l>l,m 



(21) 

Taking the Gateaux derivative of Erei[w; Q] wrt w gives the Euler-Lagrange 
equation, 

[G-2Kei]iw'^) = ^nCi;,,. (22) 

By the Fredholm Alternative [SB], equation ()22|1 either (a) has solutions 
for all values of the right hand side, or (b) has infinite number of solutions 
when the right hand side is orthogonal to the kernel of [G — 2Xrei]- All the 
eigenvalues of G are negative and form an increasing sequence, A^"^. There 
are several subcases, namely, (1) Xrei G {—oo, — |), (2) Xrei G (— ;|, oo), Xrei 
7^ ~2l(]+Ty ('^) -^^^ ~ ~ 21(1+1) ^ ^hich fall into the two broader 

classes that depend on whether or not the kernel of [G — 2Xrei] is trivial. 

Although the results in the following theorem clearly have physical sig- 
nificance, they are nonetheless easy to establish because the Euler-Lagrange 
equation is a linear equation. Figure 4 shows plots of (j^Hj) and ()24p . 

Theorem 5: 

(1) Only when Xrei £ (— oo, — |), can extremal relative vorticity in the 
form of solid- body rotation in the same direction as planetary vorticity arise 

(2) For Xrei ^ oo), Xrei 7^ ~ 21(1+1) ' cxtrcmal vorticity (if it exists), 
is solid-body rotation in the opposite direction as planetary vorticity; 

(3) Only when Xrei = ~ 21(1+1) ^ higher spherical harmonics 
ipim^ 7^ contribute to the extremal vorticity. 

The proof is divided into two natural subcases, namely, (a) when the 
kernel of [G — 2Xrei] is empty, and (b) when it is not. 

(a i) Xrei e (-CX), -i) 

In this regime, the kernel of [G — 2Xrei] is also empty, which implies the 
results: 

1 QC 

for Xrei G (-00, - -), W° = --y^ ^^^^^^O; (23) 

2 (^2 + ^Xrel J 



and. 



= kipiQ with /c ^ 0^ as A 



rel 
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W° = klpi(^ with — > CXD as Xrel — ^• 

Thus, there is a pole-hke singularity at Xrei = —1/4. This proves part (1) of 
the theorem. 

(a ii) Xrel e ( -^,0) and Xrei ^ ^iHTn) 
In this case, the following holds: 

1 1 

for Xrel e (-7,0), Xrel 7^ 



2/(/ + l)' 

(a iii) Xrel > 

Since all eigenvalues of G are negative, the kernel of [G — 2Xrei] is empty 
for Xrel > 0. Thus, the solution of is easily found by setting w° = kipiQ : 

QC 



(1 + 4:Xrel) 

This implies the result: 

= kip^Q with A; — > 0~ as A,. — > oo; 
w° —QCipiQ as A^ez — ^ 0^. 

The special solution w° = —QCtpiQ holds when the constraint Vrei is in- 
operative, and the variational problem is the unconstrained optimization of 
the energy H. Part (2) of the theorem is proved in a(ii) and a(iii). 

(b) Xrel e [-|, 0) and Xrei = -2i(TTT) 

At the bifurcation values of Xrei in this regime, the following result holds, 
as is easily ascertained: 

forA^eZ = 777^ G i--,0), / = 2,3,..., 

2Z(/ + 1) ^ 4' ^' ' ' ' 

nc ' 

W° = TT^lO^ ^ (^Im^lm- (25) 



,2 1(1+1)) ™=-' 

This proves part (3) of the theorem. QED. 
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5 Existence and properties of extremals 



The final steps of tlie Euler-Lagrange Multiplier Method compute the value 
of the multiplier Xrei by applying the fixed value of the enstrophy Qrei > 0. 
A consequence of the completion of this circle of calculations is the determi- 
nation of the physical properties of the extremal relative vorticity w'^ on the 
basis of the relative enstrophy Qrei and the spin rate fl. Figure 3 shows a 
plot of (I2ni) and (EH). 

Lemma 6: The Lagrange Multipliers Xrei of the extremals of the varia- 
tional problem Ull}) . Hl^) are given in terms of the spin rate Q > and the 
relative enstrophy Qrei > by 



A 



rel 



X. 



rel 



V Qrei 

nc 



rel 



(26) 
(27) 



Proof: Substituting the solution of the Euler-Lagrange equation into 



the constraint in Vr. 



\w 



Ol |2 



Qrei > 0, 



4(1 + 2A,,e 

gives the results for A^^^. The remaining special solutions correspond to 
a countable set of bifurcation values. QED. 

The proof of the following statement is left to the reader. 

Lemma 7: (1) The first branch of solutions in Lemma 6, 



Qreli' 



10 



(28) 



corresponding to 



\el < 



1 

■4' 



are associated with solid-body rotation in the direction of spin Vt. In terms 
of the original kinetic energy 



H 



Max 



lOJ 



^ Qrei 



;i + 8A+,) 



\ + '^Kel) 
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-r2C y Qrei- 



(29) 
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(2) The second branch of solutions 



mm 



,0 




10 



(30) 



associated with 



Kel 7^ 



'rel 



G 




/ = 2,3,..., 



2/(/ + l) 



correspond to solid-body rotation opposite the direction of spin. In this case, 
the original kinetic energy 



for given relative enstrophy Qrei md spin Vt. 

Because the Euler-Lagrange method gives only necessary conditions, to 
get sufficient conditions for the existence of constrained extremals, it is tradi- 
tional to use the so-called Direct Method of The Calculus of Variations [21] , 
|2Uj . [221, where it is shown in two parts, that (i) the unconstrained extremals 
of an augmented objective functional exist, and (ii) these unconstrained ex- 
tremals are the constrained extremals of the original objective functional. 
This classical approach is presented in a sequel The approach in this 
paper of using the infinite dimensional geometry of the energy and relative 
enstrophy manifolds, can be viewed as an alternative rigorous method for 
proving the existence of constrained extremals in the coupled fluid-sphere 
model. 

The Euler-Lagrange method gives necessary conditions for extremals of 
a constrained problem. Sufficient conditions for the existence of extremals 
can nonetheless be found from the geometrical basis of the method itself. 
An extremal must be contained in the set of points p with common tangent 
spaces of the level surfaces of H and F, but only those points p where one 
level surface remains on the same side of the other level surface in a full 
neighborhood of p correspond to constrained energy maximizers or minimiz- 
ers. Of those points, if moreover, the level surface of H is on the outside 




(31) 
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of the surface of F in a full neighborhood of PMax-, with respect to the point 
Po = {—^C, 0, 0) in the subspace of -^2(5'^) defined by (fT^. then puax is a 
constrained energy maximum. If on the other hand, the level surface of H is 
on the inside of the surface of F in a full neighborhood of Pmin, with respect 
to the point Po = {—flC, 0, 0), then pmin is a constrained energy minimum. 

All of the above Lemmas and Theorems form that part of variational 
analysis known as Necessary Conditions. We will use these facts and lemma 
3 in the proof of the following existence result which gives Sufficient Con- 
ditions for constrained extremals in terms of the geometry of the objective 
and constraint functionals. Figure 5 depicts case (1) in Theorem 8; figure 6 
depicts case (2C) and figures 7a and 7b show the retrograde solid-body state 
as constrained local minimum and maximum respectively in both case (2A) 
and (2B). 

Theorem 8: (1) The first branch of solutions w^^^ = ^/Qreii'io 
Lemmas 6 and 7 are global energy maximizers for any relative enstrophy 
Qrei and spin Vt. 

(2) For the second branch of solutions Wmin = —y/Qrei^w in Lemmas 6 
and 1, the following statements hold: 

(A) If relative enstrophy is large compared to spin, i.e., 

Qrel > 4fi^C2 (32) 

then wJJ^in is a special saddle point, in the sense that, it is a local energy 
maxima in all eig en- directions except for span{tp-^^ , where it is a local 
minima. 

(B) // relative enstrophy satisfies 

Q^C^ < Qrel < 4fi^C2 (33) 

then wjjjin is a constrained energy saddle point. 

(C) If relative enstrophy is small compared to spin, i.e., 

Qrel < ^'C^ (34) 

then wjjjj^ is a constrained energy minima. 

Proof: (1) From the eigenvalues Xim = —/(/ + !) and the fact that the 
spherical harmonics ipi^ diagonalizes the energy H, and Lemma 3 (cf. eqn. 
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((TTj)), we deduce that H is positive-definite in L2(S'^), and its level surfaces 
are infinite dimensional unbounded ellipsoids centered at Po = {—^C, 0, ...0), 



with the shortest semi-major axes (of equal lengths) 



(35) 



in span {Vio) ^11) V"!,-!} ; 



and 



all remaining semi-major axes (associated with 



(36) 



azimuthal wavenumber / greater than 1) 
have lengths L(/) that are quadratic in / 
and independent of the wavenumber m. 



The level surfaces of relative enstrophy are non-compact concentric spheres 
centered at in L^(S'^), 



Together this implies that the level surface of H lies on the outside (wrt Po) 
of the relative enstrophy level surface for fixed Qrei in a neighborhood of 
the point w%^^{Qrei) = +VQreii'io- This proves that w%j^^{Qrei) is a global 
constrained energy maximizer. 

(2) Using the same nomenclature in (D), we determine that when Qrei < 
n^C^, at the common point tf^j^ = —^Qreiipw, the relative enstrophy hyper- 
sphere lies on the outside of the energy ellipsoid wrt po] indeed both surfaces 
are locally convex with respect to their respective centers. This implies that 
wj^jj^ is a constrained energy minimum in the case ()34|) . 

When Qrei > 4i7^C^, it is natural to separate the common tangent space 
of the ellipsoid and the sphere at wj^in i^^o orthogonal parts, namely, 



/^From the property (jHSj) and ()32j). we deduce that the semi- major axes in 
span lipiQ, V'li; "^i.-i f of the energy ellipsoid at w^jn, have equal lengths 



2 _ 
2 ~ 



'^^Im = Qrei > 0. 



Im 



(a) span{V'i,±i}, (6) span / > 1}. 



(37) 



L = \- ^JqZi + ^c\ > nc. 



The enstrophy sphere at Wmin has radii of equal lengths 
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From the fact that the center of this elhpsoid is at Po = —VlCiIjiq while the 
sphere is at the origin, it follows that 

L < L^e/- 

This means that in (a) span{ipi .^i\ of the common tangent space (jHTj) at 
wjjjin, this ellipsoid is inside the sphere wrt po, which implies that w^:^^ is 
a local energy minimum there. And using ()36|) . we deduce that, in the 
orthogonal complement, part (b) of this common tangent space flHTjl at 
the ellipsoid is outside the sphere wrt po, which implies that if^in is a local 
energy maximum there. Thus, w^-^^ is a special saddle point in the case (jS21), 
in the sense that, apart from span^ipi j.^}, it is a local maximum. 

When VPC"^ < Qrei < 4i7^C^, it follows from property and the 
unboundedness property of the energy ellipsoid, that there is a critical 
value of the azimuthal wavenumber Icrit, such that in part 

(a) span{ip^^, I < Icru} 

of the common tangent space at w^J^j^, the ellipsoid is inside the sphere wrt 
Po', and in the orthogonal complement 

(6) span{i;i^, I > Icru}, 

the ellipsoid is outside the sphere wrt Po- Thus, w^jn is a constrained energy 
saddle point in the case (jH^ . QED 

6 Nonlinear Stability Analysis 

The linear stability of the steady-states WMax and 

^min CcLll be read off the 

second variation of the augmented energy functional Ej-ei- We are interested 
in proving a stronger result, namely the nonlinear or Lyapunov stability of 
the energy extremals. In principle this could be done using the geometric 
(convexity) properties of the relative enstrophy and energy manifolds, recall- 
ing that the former is a sphere in L^(5'^) and the latter is an "ellipsoid" in 
the same Hilbert space. Due to the infinite-dimensionality of L'^{S'^), and 
the consequent non-compactness of the enstrophy manifold, this is a delicate 
undertaking. 
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Instead, we propose to follow the method discussed in detail in jT7j, [TT] . 
[T^ . The first step consists of the identification of two quadratic functionals 
Qi and Q2 in the perturbation Aw such that 

Ql{Aw) < H[wMax + Aw] - H[wMax] - DHMax " Aw 
Q2{Aw) < T[wMax + Aw]-T[wMax]- D'^Max- 

Let {aim} denote the perturbation Fourier coefficient based on the steady- 
state WMax = VQreii^ 10, that is, the steady state is given by 

aio{wMax) = \J Qreh 

(^imiwMax) = for all other modes. 
It is natural and convenient to choose the positive definite quadratic forms 
Qi({a/m}) 



Q2{,{aim}) 

Then, Qi + Q2 is positive definite, and is clearly a norm for the perturbation 
Aw since Q2{Aw) = ||Aw||2. Indeed, it is the energy-enstrophy norm of 
Arnold. 

From the fact that at steady-state WMax, D{H + T){wMax) ■ Aw = 0, it 
follows that 

(Ql + Q2){{aimit)}) = iH + T)[wMax + Aw{t)]-{H + T)[wMax] 

= (H + T) [wMax + Au;(0)] -(H + T) [wMax] 

< AE{0), 

that is, the sum Qi + Q2 is conserved along trajectories of the BVE. Thus, 
we conclude that trajectories that start near WMax will remain nearby for all 
time, that is, WMax has been shown to be Lyapunov stable. 

The same arguments go through for the energy minimizer Wmin when 
the relative enstrophy Q^ei < ^"^C^- And for the saddle point w^i^ when 
Qrei > 4n^C^, a modified argument where perturbations are restricted to 



= H[WMax + Aw] - H[wMax] " DHMax " Aw 

A 9 \ — ^' 

^ ^ Zm^lO 

= T[wMax + Aw] -T[wMax] - DTMax ■ Aw 

= E«L>o. 
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be in the orthogonal complement of span{ilJiQ j.^} in L^(S'^), proves that 
Wmin has some degree of nonlinear stability. This completes the proof of 

Theorem 9: The global energy maximizer WMax is Lyapunov stable un- 
der all conditions of spin, energy and relative enstrophy. The constrained 
energy minimizer Wmin is Lyapunov stable in the low energy, low enstrophy 
regime, that is, when Qrei < fl^C^- If Qrei > 4:^^C^ (high enstrophy, high 
energy regime), then the saddle point uVin is nonlinearly stable to all per- 
turbations with azim,uthal wavenumber I > 2. In the intermediate regime, 
Q'^C'^ < Qrei < 4:0,"^ C^, Wiain is Unstable. 

7 Discussion and conclusions 

Prom the point of view of applications to planetary atmospheres, the physical 
consequences of theorems 5, 8 and 9 is embodied in the following over-arching 

statement: 

There is a relative enstrophy threshold value Vt^C^such that below it, the 
constrained energy extremals consist of the counter-rotating local energy min- 
imizer state = —y/Qreii^ 10 c^HC^ the pro-rotating global energy maximizer 
state w^^^ — -t-VQre/V'io- When relative enstrophy exceeds this value, the 
pro-rotating global energy maximizer state tf^^^ = +^/Qrei'4'io Persists but 
the counter-rotating state w^^^^^ becomes a constrained saddle point. 

In other words, the counter-rotating state is a nonlinearly stable energy 
extremal only for low relative enstrophy and low kinetic energy. The nonlin- 
early stable counter-rotating steady state w^^^ — —y/Qreii'io changes stability 
to a constrained saddle point via a saddle-node bifurcation when the rate of 
spin fl decreases past the threshold value 

Q _ VQrel 

°~ C ' 

with relative enstrophy fixed at Qrei- When the relative enstropy is large 
enough compared to the spin rate Q > 0, that is, Qrei > il^C^, the only 
nonlinearly stable steady state in the model is the global energy maximizer 
'^Max — VQrei'ipiOi with the cavcat that for Qrei > 4^2^(7^, the counter- 
rotating state Wmin i^ nonlinearly stable to zonally- symmetric perturbations. 

Thus, for any relative enstrophy and any given spin rate Q > 0, a super- 
rotational Lyapunov stable steady state exists in the coupled model. When 
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Qrei < ^^C^, both the global energy maximizer and global energy minimizer 
are Lyapunov stable solid-body rotational states. 

The range of phenomena associated with the study of atmospheric super- 
rotation such as pertaining to the Venusian atmosphere (cf. [13^ and refer- 
ences therein) is very complex. While we do not claim that our theory and 
results explain the whole complex range of phenomena associated with super- 
rotation, by going to a model of a barotropic fluid coupled to an infinitely 
massive rotating sphere - the fluid component exchanges angular momentum 
as well as kinetic energy with the sphere - we have obtained rigorous results 
that are directly related to atmospheric super-rotation in both the qualitative 
and quantitative sense. 

We will compare our results here with those from simulations of the Gen- 
eral Circulation Model (GCM) as well as those obtained from intermediate 
level models. The starting point of the body of work based on the GCM is a 
multi-layer stratified rotating flow model. Furthermore, the system modelled 
is usually a damped driven one with differential solar radiation and Ekman 
damping. The best phenomenological results to date - using the GCM - are 
therefore necessarily based on intensive numerical simulations. 

For the purpose of this comparison, numerical simulations of the ter- 
restrial GCM show that super (sub)-rotational flows behave like saturated 
asymptotic steady states of a complex damped driven system (cf. UH] and 
ref. therein). The super- rotating state appears to be associated only with 
sufficiently energetic flows in these simulations. Moreover, the basin of at- 
traction of these flows can be reasonably characterized by only a few initial 
quantities such as energy and angular momentum, and the super (sub) - 
rotational end states appear to be quite robust in a numerical sense. The 
super-rotational state is taken to be the one where the fluid component has 
positive net angular momentum - in the direction of the sphere's spin - rel- 
ative to the frame in which the sphere is fixed. The vertically averaged 
barotropic component of the flow in a damped driven multi-layered atmo- 
spheric system therefore have asymptotic quasi-steady states which are super 
and sub-rotational states near to the basic zonal states ± ipio [131 • R-ecent 
numerics (cf. ^3] and references therein) confirm the super-rotational state 
is more relevant to the atmosphere of Venus and is associated only with suf- 
ficiently high kinetic energy; the sub-rotational state arises only for lower 
energies in for instance, the atmosphere of Titan 13]. 

A more direct comparison can be made between our detailed results and 
those reported by Yoden and Yamada 0. In particular, they found robust 
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relaxation to a prograde solid-body state for small to intermediate values of 
the planetary spin, and to a retrograde solid-body state when the planetary 
spin is large relative to the kinetic energy in the barotropic flow. The first is 
consistent with our first prediction that the prograde sohd-body state arises 
only when the rest frame kinetic energy of the flow is high enough relative 
to the planetary spin, or equivalently, for a given range of kinetic energy, it 
is allowed only for relatively small planetary spins. The second is partially 
consistent with our second prediction that the retrograde solid-body state is 
nonlinearly stable only for planetary spins that are large in comparison to the 
relative enstrophy, if one allows for the fact that the large anticyclonic polar 
vortex state is a superposition of mainly the retrograde solid-body state and 
some zonally symmetric spherical harmonics ipiQ with wavenumber / < 5. 

A related point vortex formulation for the inviscid dynamics of a rotating 
barotropic vorticity field has been recently reported by Newton et al in 
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Figure 1: Graph of energy H vs. coordinates k 



a 10 of extremals as in 



fig. 1 




Figure 2: Energy H - relative enstrophy Qrei space; black region denotes un- 
permitted values; gray region denotes non-extremal permitted values; curves 
denote values at extremals as in (fT^. 
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Figure 4: Graph of extremal coordinates k vs. Lagrange Multipliers X^^i for 
fixed spin > as in (j22D and (1^ . 
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Figure 5: Projections of energy ellipsoid and enstrophy sphere showing 
common tangent at global maximizer w^^^ when energy exceeds H^. 

fig. 5 
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Figure 6: Projections of energy ellipsoid and enstrophy sphere showing 
common tangent at local minimizer w^^^ when fl34j) holds. 
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Figure 7: (a) Projections of energy ellipsoid and enstrophy sphere showing 
the saddle point w'^^^ as local minimum when or holds; (b) Projec- 
tions of energy ellipsoid and enstrophy sphere showing the saddle point w'^^^ 
as local maximum 

fig. 7(a) 
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fig. 7(b) 
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List of Figure Captions: 

Figure 1. Graph of energy H vs. coordinates k = aio of extremals as in 
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Figure 2. Energy H - relative enstrophy Qrei space; black region de- 
notes unpermitted values; gray region denotes non-extremal permitted val- 
ues; curves denote values at extremals as in (fT^ . 

Figure 3. Graph of Lagrange Multipliers A^^^ vs. square- root of relative 
enstrophy, Qrei for fixed spin > as in and (P7j). 

Figure 4. Graph of extremal coordinates k vs. Lagrange Multipliers Xf^i 
for fixed spin > as in and 

Figure 5. Projections of energy ellipsoid and enstrophy sphere showing 
the common tangent at global maximizer when energy exceeds He- 

Figure 6. Projections of energy ellipsoid and enstrophy sphere showing 
the common tangent at local minimizer w^j^ when fl34p holds. 

Figure 7a. Projections of energy ellipsoid and enstrophy sphere showing 
the saddle point w^J^i^ as local minimum when (jHHj) or ((221) holds. 

Figure 7b. Projections of energy ellipsoid and enstrophy sphere showing 
the saddle point w'^^^ as local maximum. 
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